Emergence of quantum chaos in finite interacting Fermi systems 
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We study the level spacing statistics P(s) in many-body Fermi systems and determine a critical 
two-body interaction strength U c at which a crossover from Poisson to Wigner-Dyson statistics takes 
place. Near the Fermi level the results allow to find a critical temperature T c h above which quantum 
chaos and thermalization set in. 
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The Random Matrix Theory (RMT) was developed to 
explain the general properties of complex energy spec- 
tra in many-body interacting systems such as heavy nu- 
clei, many electron atoms and molecules Later, it 
has found many other successful applications in different 
physical systems. Among the most recent of them we 
can quote models of quantum chaos where RMT appears 
due to the classically chaotic but deterministic under- 
lying dynamics 0]. One of the most direct indications 
of the emergence of quantum chaos is the transition of 
the level spacing statistics P(s) from Poisson to Wigner- 
Dyson (WD) distribution. This property has been widely 
used to detect the transition from integrability to chaos 
not only in systems with few degrees of freedom 0j but 
also in solid-state models with many interacting electrons 
||. It was also applied to determine the derealization 
threshold in noninteracting disordered systems [Q . 

While the conditions for the appearance of the WD 
distribution in noninteracting systems is qualitatively 
well understood the situation is more intricate in pres- 
ence of interaction. Indeed, in this case the size of the 
total Hamiltonian matrix grows exponentially with the 
number of particles and it becomes very sparse as a re- 
sult of the two-body nature of the interaction. Due to 
that it was initially not obvious whether switching on 
the interaction would lead to the WD statistics. To 
study this problem a two-body random interaction model 
(TBRIM) had been proposed This model consists 

of n fermions which can occupy m unperturbed energy or- 
bitals with mean one-particle level spacing A. The mul- 
tiparticle states are coupled by two-body random transi- 
tion matrix elements of typical strength U. It was found 
that a sufficiently strong U leads to a level mixing and 
appearance of WD statistics. Very recently the interest 
for this model has been renewed and its statistical prop- 
erties were investigated in more details JtJ . This raise of 
interest was stimulated by the understanding that many 
statistical properties of real physical systems such as the 
rare-earth Ce atom M and the 28 Si nucleus PJlC|] are 
well described by the TBRIM. In addition this model is 
quite similar to the s-d shell model used for a description 
of complex nuclei P,fL0[ . Since interaction is generically 



of two-body nature it is reasonable to assume that this 
model will be also useful for a description of interacting 
electrons in clusters [ pT[ and mesoscopic quantum dots 

& 

While the statistical properties of the TBRIM were 
studied in some details, surprisingly, the most important 
question of the critical interaction strength U c at which 
the WD level spacing statistics sets in was omitted. Ap- 
parently the reason for this is based on the common lore 
in nuclear physics that the level density grows exponen- 
tially with the number of particles and therefore an ex- 
ponentially small interaction is sufficient to mix nearby 
levels gjlO . However recent estimates on few-particle 
models (n = 2,3,4) showed that in spite of the high 
many-body density of states, only an interaction strength 
comparable to the two-particle level spacings can give a 
level mixing jl3|,[l4| . Therefore the dependence of U c on 
the number of particles and orbitals as well as the excita- 
tion energy should still be determined. This is the main 
purpose of this paper. The above border in U is physi- 
cally very important. Indeed for U < U c levels are not 
mixed by interaction and hence the system is not ther- 
malized. Consequently the occupation numbers are not 
described by the Fermi-Dirac statistics. On the contrary 
a sufficiently strong interaction leads to thermalization 
as it has been seen in numerical simulations J7p|jl0|] . 

To study the effect of interaction on the spectral prop- 
erties of finite Fermi systems we used the TBRI model 
described in M. It consists of n particles distributed 
over to orbitals with energies e m i , m' = 1,2, ...m. 
These energies are randomly distributed over the in- 
terval [0, 777.] with average spacing A = 1. The total 
number of multiparticle states is N = m!/(n!(m — 77)!). 
They are coupled by random two-body transition ma- 
trix elements distributed in the interval [—U,U]. Due 
to the two-body nature of the interaction, only states 
differing by at most two one-particle indices are cou- 
pled. As a result each multiparticle state is coupled with 
K = 1 + n(m — 77) + 77(77 — l)(m — 77) (777 — 77 — l)/4 states 
. All these transitions occur inside a two-body energy 
interval B — 2m — 4 around the energy of an initial mul- 
tiparticle state. For large m and 77, the number of tran- 
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sitions K is much smaller than the size of the matrix N 
but is much larger than the number of different two-body 
matrix elements N2 ~ m 2 /2. The total energy of the sys- 
tem varies from the ground-state value E g w n 2 A/2 to 
the maximal value E t w mnA — E g and the Fermi energy 
is ep ~ nA. The typical level spacing in the middle of 
the spectrum at E h « (E t +E g )/2 is A„ sa (E t -E g )/N. 
Let us first discuss the situation at high energies E ~ 
where all if transitions are energetically allowed. In 
this case the density of directly coupled states is p c ~ 
K/B because all transitions take place inside the two- 
body energy band B. According to perturbation theory 
these levels will be mixed when the transition matrix 
element U between them becomes of the order of the 
corresponding spacing A c = l/p c . This determines the 
critical coupling U c 
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Here, we introduced the two-particle density p 2 ~ 
N2/B « m/4 assuming m 3> n 3> 1 and a numerical 
constant C to be determined. For U <C U c the pertur- 
bation theory works, levels are not mixed and P{s) is 
close to the Poisson distribution. For U > U c we expect 
a strong mixing of levels not only on a scale A c but on a 
much smaller scale A„ . There are few arguments in favor 
of this statement. The first of them is based on the re- 
sults for few-particle systems (n — 2, 3, 4) (l3). According 
to jD^, the effective transition matrix element between 
nearby levels in high orders of perturbation theory be- 
comes comparable to A„ when the first-order transition 
mixes directly coupled states (U > U c ). Recently the 
same conclusion was drawn in |l4| . The second argument 
is based on an analogy with superimposed band random 
matrices (SBRM) with strongly fluctuating diagonal el- 
ements 15 - 18 1 . There it was shown that for sufficiently 
large band (number of nonzero diagonals 2b + 1 3> y/N) 
the eigenstates are extended over the whole matrix size 
N and P(s) has the WD form if the transition matrix 
elements are larger than the energy spacing between di- 
rectly coupled states. This condition is rather similar to 
the above border (jl|). 

To check the prediction (Q), we numerically computed 
P(s) in the middle of the spectrum of the TBRIM (keep- 
ing only ±25% of the levels around Eh) for n < 8 and 
m < 80 at various interaction strengths U. Up to 5000 
different realizations of disorder have been used to ob- 
tain the total spacing statistics N s w 30000. A typical 
example of the transition from Poisson to WD statistics 
is shown in Fig. 1. As expected the level repulsion dis- 
appears at small U while for large U the distribution 
approaches the WD form. To characterize this transi- 
tion we computed for each distribution P(s) the value 
V = C(P(s) - PwD(s))ds/f* (P P (s) - P WD (s))ds. 
Here Pp(s) and Pwd{s) are the Poisson and the WD 



distributions respectively and sq — 0.4729... is their inter- 
section point. In this way rj varies from 1 (P(s) = Pp(s)) 
to (P{s) = P\vd(s))- We determined the critical inter- 
action strength U c by the condition rj(U c ) = r) c = 0.3. 
The choice of r\ c influences only the numerical factor C 
in (|l]). We note that this rj- value is close to the value 
Va = 0.215 corresponding to P(s) at the Anderson tran- 
sition in 3d (in a criterion slightly different from 
ours was used). 
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FIG. 1. Transition from Poisson to Wigner-Dyson statistics 
in the TBRIM for m = 12, n = 6 : U/A = 0.01 and r) = 0.93 
(+); U/A = 0.055 and 77 = 0.3 (•); U/A = 0.13 and 77 = 0.063 
(x). Full lines show the Poisson and the Wigner-Dyson dis- 
tributions. Insert shows P(s) at fixed rj — 0.3 for half-filling 
v = n/m = 0.5 and n = 4 (dotted line), n = 5 (dashed line), 
n = 6 (long dashed line) and n = 7 (dotted-dashed line) . 

The fact that the concrete choice of r\ c is not crucial 
is also confirmed by Fig. 2 which shows the existence of 
a scaling rj = rj(U/U c ). Indeed the numerical data in a 
large parameter range demonstrate the existence of one 
scaling curve (Fig. 2). This scaling is very similar to the 
one observed in the SBRM models p5|— p~8[ . It also clearly 
shows that the situation in our model is qualitatively 
different from the ^-scaling in the solid-state models with 
Anderson transition. There, in the limit of large system 
size, only three values r/ = 1 (localized phase), r\ = 
(delocalized) and 77 = r/A (at the transition) are possible 
]4j. On the contrary in our case the scaling function 
varies smoothly from 1 to with the rescaled transition 
matrix element U/U c for different system sizes N which 
varied over more than two orders of magnitude. We relate 
this qualitative difference between the two models to the 
fact that in the TBRIM all orbitals are coupled by direct 
transitions whereas in the Anderson model, the hopping 
couples only nearby sites. Due to that the TBRIM is 
more similar to the SBRM models with broad band where 
many states are directly coupled. 
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FIG. 2. Dependence of 77 on the rescaled interaction 
strength U/U c for 2 < n < 8, 4 < m < 80, 1/40 < z/ < 1/2 
and 0.02 < U c < 0.2 (diamonds). Open circles show the scal- 
ing close to the Fermi level (see text). The straight line marks 
n = -q c = 0.3. 

The condition for the critical U c (n c — 0.3) allows to 
check the theoretical prediction (|l]) . The numerical data 
for which the number of direct transitions varies over 
more than two orders of magnitude are presented in Fig. 
3. They give a clear confirmation of the estimate ([!]) 
giving C w 0.58. The results of Figs. 1-3 show that for 
U > U c from (|l|) all nearby levels are mixed by two-body 
interaction and P(s) converges to the RMT result with 
WD distribution. We stress that for large m and n, the 
value of U c remains parametrically much larger than the 
multiparticle spacing A„. 
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FIG. 3. Dependence of the rescaled critical interaction 
strength U c /B, above which P(s) becomes close to the 
Wigner-Dyson statistics, on the number of directly coupled 
states K for 4 < m < 80 and 1/40 < v < 1/2. The line shows 
the theory (1) with C = 0.58. 



So far the results were obtained in the middle of the 
energy spectrum Eh where all K direct transitions are 
energetically allowed and effectively work. The situation 
becomes quite different close to the Fermi level. There, 
the estimate (|l|) should be modified in the following way. 
First we should take into account that the density of effec- 
tively coupled two-particle states p2ef becomes energy- 
dependent so that y02e/(e) ~ e/A 2 [19 2(J. Secondly the 
number of effectively interacting particles is also changed 
close to the Fermi level. Indeed as it is well known, at 
a temperature T, only Sn ~ Tnjtp ~ T/A > 1 parti- 
cles interact near the Fermi surface. At this excitation 
energy e ~ T < ejr, the density of two-particle states is 
P2ef ~ T/A 2 . By replacing in ([!]) n by Sn and p2 by 
P2ef we obtain that at a given interaction strength the 
levels become mixed and P(s) takes the WD form at a 
temperature higher than the critical T c h given by 
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where C\ is a numerical constant. The conditions of 
validity of this equation are T c h > A (Sn > 1) and 
T c h < cf = riA that corresponds to n~ 3 < U/A < 1. It 
is also assumed that the WD statistics implies thermal- 
ization with Fermi-Dirac statistics. Such a conjecture 
looks quite natural, since the quantum chaos should be 
related with excitation of many unperturbed modes and 
mixing. Also without mixing of nearby levels and WD 
statistics the thermalization is not possible since gener- 
ally the Poisson distribution indicates an existence of un- 
coupled parts in the whole system. As a result the ther- 
malization does not exist below T c h- 

Since near the Fermi level the total system energy 
counted from E g is 8E = E — E g w TSn, the relation 
(g) implies that the thermalization takes place only for 
eigenstates with eigenenergies E\ = E g + SE so that 



SE > SE ch ~ C 2 A(A/U) 
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The above restriction for U requires 1 < SE/A < n 2 . 
This result shows that the ^-parameter should depend 
on the excitation energy. Indeed, our numerical data, 
extracted from P(s) computed in a small energy inter- 
val near a fixed SE, clearly show that n decreases with 
increasing excitation energy SE (Fig. 4). Using the 
relation (^|) we can determine for a given SE an effec- 
tive U c value being U c = CfA(A/SE) 3 / 2 . The con- 
dition rj(SE) = rjc = 0.3 for the data of Fig. 4 at 
n = 6,m = 12, U/A = 0.147 gives d » 1.08. With 
the value C\ = 1.08 and the above dependence of U c 
on SE we can check if the data of Fig. 4 will follow 
the general scaling law of Fig. 2. For that in Fig. 2 
we plot the ^-values of Fig. 4 vs. the ratio U/U c with 
U c = 1.26((5£)- 3 / 2 , Ci = 1.08 and A = 1 (open circles). 
The fact that these data follow the scaling curve confirms 
the theoretical estimates (||), (||) for the thermalization 
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border. The direct check of the dependence of SE c h on 
U (insert in Fig. 4) also confirms the prediction (3). 
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FIG. 4. Dependence of r\ on the rescaled excitation energy 
5E/A for n = 6,m = 12 and U/A = 0.147 (o). The straight 
line marks rj = r\ c = 0.3. Insert gives the numerically found 
dependence of SE c h with rj = r\ c = 0.3 on U (diamonds), the 
straight line shows the theory (3) with C\ = 1.08. 

The obtained estimates for the quantum chaos bor- 
der (||), (||) can be applied to different finite interacting 
Fermi systems such as complex nuclei with residual inter- 
action, atoms and molecules, clusters and quantum dots. 
Here we briefly discuss the case of metallic quantum dots 
JTH . In this case the interparticle interaction is relatively 
weak so that U/A ~ 1/g with g = E c /A 3> 1 being the 
conductance of the dot and E c the Thouless energy [glj . 
According to (|^) the thermalization will take place above 
the excitation energy 6E ch ~ Ag 2 / 3 . This is in a satisfac- 
tory agreement with the experimental results [ [l2[ where 
a dense spectrum of excitations in dots with g ~ 100 ap- 
pears at excitation energies 8E c h ~ 10A. We note that 
our border for thermalization and chaos 5E c h is higher 
than the border for quasiparticle disintegration on many 
modes SEp ~ Ag 1 / 2 proposed in p2| , p3[ . In our opinion 
the parametrically different dependence on g suggested 
in |2^,|2j| appears because the effect of energy redistri- 
bution between many excited modes was neglected while 
the derivation of estimates (||), (g|) shows that it plays 
an important role. In addition, in the relations similar 
to (||) the authors of |22|,^3| in fact used the first power 
of n, instead of n 2 , that according to our numerical data 
(Fig. 3) does not correspond to the regime with many 
excited modes. 

In conclusion, our numerical results and analytical es- 
timates allowed to determine the border for emergence 
of quantum chaos and thermalization in finite interact- 
ing Fermi systems. Further investigation of this crossover 
in real systems is highly desirable. 
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